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ABSTRACT 

We report the detection of 21-cni and molecular hydrogen absorption lines in the 
same damped Lyman-a system (with log N{I1 i)=:21.36±0.10) at Zabs^S. 17447 to- 
wards SDSS J133724.69+315254.55 (zcm ~3.174). We estimate the spin temperature 
of the gas to be, Ts = 600l^gg K, intermediate between the expected values for cold 
and warm neutral media. This suggests that the H i absorption originates from a 
mixture of different phases. The total molecular fraction is low, /h2 = 10~^, and H2 
rotational level populations are not in equilibrium. The average abundance of the a- 
elements is, [S/H]=— 1.45 ± 0.22. Nitrogen and iron are found underabundant with 
respect to a-elements by '^1.0 dex and ^0.5 dex respectively. Using photoionization 
models we conclude that the gas, of mean density, nn ^2 cm~'^, is located more than 
270 kpc away from the QSO. While the position of 21-cm absorption line coincides 
with the H2 velocity profile, its centroid is shifted by ~2.7±1.0 km s^^ with respect 
to the redshift measured from the H2 lines. However, the position of the strongest 
metal absorption component matches the position of the 21-cm absorption line within 
0.5 km s~^. From this, we constrain the variation of the combination of fundamental 
constants x = a'^Gp/fi, Ax/x — — (1.7±1.7) x 10~^. This system is unique as we can at 
the same time have an independent constrain on fi using H2 lines. However, as the H2 
column density is low, only Werner band absorption lines are seen and, unfortunately, 
the range of sensitivity coefficients is too narrow to provide a stringent constraint: 
AyLt//i < 4.0 X 10"^. The Uhraviolet and Visual Echelle Spectrograph (UVES) spec- 
trum reveals another DLA at Zabs= 3.16768 with log A^(H i) = 20.41±0.15 and low 
metallicity, [Si/H] = -2.68 ± 0.11, in which [0/C] - 0.18 ± 0.18 and [0/Si] - 0. 
This shows that even in the very early stages of chemical evolution, the carbon or 
silicon to oxygen ratios can be close to solar. Using Voigt profile fitting we derive 
log(A^(D i)/A^(H i)) = -(4.93 ± 0.15) in this system. This is a factor of two smaller 
than the value expected from the best fitted value of fib from the Wilkinson Microwave 
Anisotropy Probe (WMAP) 5 year data. This confirms the presence of astration of 
deuterium even at very low metallicity. 

Key words: Galaxies: abundances - Quasars: absorption lines - Quasars: individual: 
SDSS J133724.69+315254.55 
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1 INTRODUCTION 

Damped Lyman-a systems (DLAs) are the highest H I col- 
umn density absorbers seen in QSO spectra, with A^(H l)> 
10^" cm~^, over cosmological time-scales. These absorbers 
trac e the bulk of neutral hydr ogen at 2 < 2: < 3 (for example 
see, iNoterdaeme et al]|2009bl ) and have long been identified 
as revealing the interstellar medium of the high- redshift pre- 
cursors of present day galaxies (for a review see, IWolfe et all 
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The physical conditions in the diffuse interstellar gas 
are influenced by a number of physical processes including 
in-situ star-formation, cosmic ray interaction, photoelectric 
heating by dust as well as mechanical energy input from both 
impulsive disturbances like supernova explosions and steady 
injection of energy in the form of stellar winds. Therefore, 
studies of the physical conditions in DLAs provide a probe of 
the cosmological evolution of normal galaxies selected solely 
based on H l cross-section. 

Our understanding of physical conditions in DLAs is 
primarily derived from optical absorption-line studies, in- 
volving the detection of low- ionization metal transitions and, 
in some cases, mole c ular abso rption [e.g. iPetitiean et al.l 
l|2000l ). lLedoux et all (|2003|) a ndlNoterdaeme et all (l2008al), 
for H2 , and ISrianand et al] (120081 ') and iNoterdaeme et al.l 
l|2009al ) for CO]. The gas in which H2 is detected, has typi- 
cal temperature and density of 150±80 K and nn = 10—200 
cm~^ respectively. These conditions correspond to those pre- 
vailing in a cold neutral medium (CNM) embedded in a ra- 
diation field of mod erate intensity origina ting from local star 
formation activity l|Srianand et al.luOOSl '). 

Detecting 21-cm absorption in DLAs is one complemen- 
tary way to probe the physical conditions in the absorbing 
gas. If detected, the 21-cm optical depth, in conjunction with 
the observed N{}i i) from op tical data, gives the spin tem - 
perature of the H I gas (see iKanekar fc C hcngaluJ l2003h . 
It is widely believed that the latter is a good tracer of the 
kinetic temperature. The width of the 21-cm line observed 
at very high spectral resolution can also yield a direct mea- 
surement (or give a stringent upper limit on) of the kinetic 
temperature (as, for example, in the 2:a bs= 1.3603 system 
towards J 2340-0053. ICupta et allbOOgl '). This can be com- 
bined with measurements of C I, C 11, Si 1 1 and /or O I fine- 
struc ture lines to derive the gas density (|Heinmiiller et al] 
I2OO6I ). 

By measuring the relative positions of absorption lines 
it is possible to constrain the time variation of several dimen- 
sionless constants. Accurately measured wavelengths of H2 
absorption lines allow one to co nstrain the variation of the 
proto n-to-electron mass ratio (jVarshalovich fc Levshakovl 
I1993I ). The position of the 21-cm absorption line can be com- 
pared to those of the corresponding metal line components 
to probe the variations of a combination, x = a^Gp/fj,, of the 
fine structure constant (a) , the proton-to-ele ctron mass ratio 
(/x) a nd the proton gyromag netic factor fGp. lTubbs fc Wolfd 
1 19801 ). Simultaneous detection of H2, 21-cm and metal ab- 
sorption line in the same system can, in principle, allow one 
to lift part of the degeneracy between the variations of these 
different constants. 

Despite obvious advantages, till recently, there was 
only one system (at ^abs ~ 1.7765 towards Q 1331-1-170, 
ICui et cil]|2005l ) known to produce both H2 and 21-cm ab- 
sorption lines. Due to the low redshift of the system, how- 
ever, it is not possible to cover the H2 Lyman and Werner 
band absorption with high resolution spectroscopy at high 
enough signal-to-noise ratio. This prevents one from per- 
forming a detailed analysis of the cloud. In the course of 
our on-going GBT/GMRT survey for 21-cm absorption in 
high-z DLAs selected from the SDSS we discovered a DLA 
system, exhibiting both 21-cm and H2 absorption lines, 
at 2abs = 3.17447 towards J133724.69-H315254.55 (called 
J 1337+3152 hereafter). This system is unique as both 21- 
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Figure 1. Normalised GMRT spectra showing the detection of 
z = 3.174480 21-cm absorption toward J 1337-1-3152. Top panel: 
low resolution spectrum (6.88 km/s/channel) observed on Jan- 
uary 13, 2009. Middle and bottom panels: high resolution spec- 
trum (1.72 km/s/channel) together with, respectively, the best 
single or two gaussian component fit models. Residuals (shifted 
by 0.6 along the y-axis) are also shown. 



cm absorption and H2 are extremely rare in DLAs at z >3. 
In addition, there is another DLA at z^hs = 3.16768 along 
the same line-of-sight (at ~500 km s~^ away from the first 
DLA) with very low metallicity. In this paper, we present a 
detailed analysis of both DLAs. 



2 DETAILS OF OBSERVATIONS 

The quasar J1337-I-3158 was observed with the GMRT on 
January 13, 2009 to search for 21-cm absorption in the 
DLA at 2:abs=3. 17448. Observations were performed in a 
bandwidth of 1 MHz split into 128 channels. We observed 
the nearby standard flux density calibrator 3C 286 every 
45 minutes for the amplitude, phase and bandpass calibra- 
tion. A total of 6.2 hrs of data were acquired on-source in 
both the circular polarization channels, RR and LL. Data 
were reduced using the NRAO Astronomical Image Pro- 
cessing Sy stem (AIPS) f o llowin g standard procedures as de- 
scribed in iGupta et al.l (|2006l ) . Special care was taken to 
exclude the baselines and timestamps affected by radio fre- 
quency interferences (RFI). An absorption feature confined 
in a single spectral resolution channel (6.88 kms~^) is de- 
tected consistently in both polarisations (see top panel of 
Fig.[T}. In the Stokes I spectrum, the absorption corresponds 
to an integrated 21-cm optical depth, J rdv = 1.55±0.19 
km s"\ The quasar was reobserved on March 17, 2009 
during 7.8 hrs (on-source time) at a higher frequency res- 
olution to confirm the 21-cm absorption. A bandwidth of 
0.25 MHz split into 128 channels, i.e. a spectral resolution 
of 1.72 kms~^, was used. With respect to the first observing 
run, this second epoch observation corresponds to an over- 
all shift in the observing frequency of 18.2 km s~^ due to 
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the heliocentric motion of the Earth. The absorption profile 
was observed to be consistent with this velocity shift. We 
used AIPS task CVEL to correct the observed data for the 
Earth's motion and rotation. The higher resolution Stokes 1 
spectrum (with rms of 2.0mJybeam~^ channeP^) from the 
second epoch observation together with the fits to the If I 
absorption are shown in the middle and bottom panels of 

Fig.m 

We also observed J 1337+3152 wi th the Ultravi o let an d 
Visual Echelle spectrograph (UVES; lOekker et aP (|20o3)) 
in visitor mode on March 27 and April 27, 2009. We used 
two settings with two different central wavelengths in the 
blue (390nm and 437nm) and red (600nm and 770nm) arms 
to cover the wavelength range 3300—9590 A with a single 
gap between 7590 and 7769 A. The CCD pixels were binned 
2x2 and the slit width adjusted to I", yielding a resolv- 
ing power of f? = 45000 under median seeing conditions of 
~ O'.'Q. The total exposure time on source is 4 hrs. The 
data were reduced using the UVES pipeline (|Ballester et al.l 
l2000l 'l. which is available in the context of the ESQ MIDAS 
data reduction system. Wavelengths were rebinned to the 
heliocentric rest frame and individual scientific exposures 
co-added using a sliding window and weighting the signal by 
the inverse of the variance in each pixel. The resultant com- 
bined spectrum has a SNR in the range 10-20 per pixel. We 
fitted the quasar continuum using low order polynomials in 
line free continuum regions. We analysed the spectrum using 
standard Voig t profile fitting procedures. The atomic data 
used are from iMortonI (|2003h unless otherwise specified. In 
the f ollowing, we adopt th e solar photospheric abundances 
from lAsplund et al.1 l|2009D . rest wave lengths an d osci llator 
strengths of H2 tra nsitions from lUbachs et all (|2007l ') and 
lAbgrall et al.1 ij 19941 ') respectively. 



3 HIGH iV(H I) SYSTEMS 

There are three high-7V(H l) absorption systems within 
2500 km s"'^ along this line of sight at Zabs ~ 3.17448 
{v = km s"^), 3.16768 (v = -488 km s"^) and 3.1408 
{v = —2420 km s~^). The corresponding Lyman series ab- 
sorption lines are shown on a velocity scale in Fig. [21 The 
Lyman-Q profile is dominated by the z^hs ~ 3.17448 system. 
As the red wing of this absorption is affected by the QSO 
Lyman-Q emission, we do not use it to derive N{}i l) but 
use instead Lyman- /3 and other higher order lines. 

The Lyman-^ line of the Zabs = 3.17448 and 3.16768 
systems (referred to as system- 1 and system-2 from now 
on) are blended together. Some of the higher order Lyman 
lines of system- 1 are blended with lines from the sub-DLA 

at Zabs= 3.1408 (-• 2420 km s"^). The main constraint 

on log A'^(II l) in system-1 comes from the red wing of the 
Lyman-/3 absorption and the core of the Lyman-7 absorp- 
tion. In the case of system-2, the blue wing of Lyman-/? and 
the red wing of Lyman-7 are used for the same purpose. We 
obtain log N{a l)= 21.36±0.10 and 20.41±0.15 for system- 
1 and system-2 respectively. Note that errors come mostly 
from continuum placement uncertainties. We discuss the fits 
to the H I lines for system-2 in greater detail below (see 
Sect. 5). In both cases, the unblended wings of the higher 
Lyman series lines are used to constrain the velocity disper- 
sion. 



In addition, based on the absence of damping wings 
in Lyman-/? and the blue wing of the Lyman-Q absorption 
lines we derive an upper limit of 10^® cm~'^ for A'^(II l) at 
Zabs= 3.1408 (Ai; ~ — 2420 km s~^). We cannot perform 
a more precise measurement as all the available transitions 
are saturated, hiding the component structure. We therefore 
do not discuss this system any further. 



4 ANALYSIS OF SYSTEM-1 (zabs= 3.17448) 

Weak II2 and strong 21-cm absorptions are seen associated 
with this system in the UVES and GMRT spectra, respec- 
tively. 

4.1 21-cm absorption 

The GMRT spectrum covers only the expected 21-cm range 
for system-1. In the higher resolution GMRT spectrum the 
absorption is spread over 7-8 channels and can be adequately 
fitted with a single Gaussian component of peak optical 
depth r = 0.40±0.02 located at 340.2593±0.0002 MHz (or 
Zabs = 3.174480±0.000002) and FWHM = 5.08±0.34kms-^ 
(see middle panel in Fig. [T]). 

The gaussian fit was carried out using an IDL code 
based on MPFIT, which performs y^-minimisa tion by 
Levenberg-Marquardt techniques (|Markwardtl2009l 'l . We use 
the rms flux measured in the line free channels as typical er- 
ror in each channel. The best fit has a reduced = 1.02. 
The FWHM of the fitted Gaussian component corresponds 
to a kinetic temperature of Tk < 564±76 K. 

Motivated by the fact that the centroid of the H2 ab- 
sorption does not coincide with the centroid of the 21- 
cm absorption (see below), we performed a fit with two 
components. The result of this fit is shown in the bot- 
tom panel of Fig [T] The best fitted two components have 
2abs = 3.174492±0.000007 and 3. 174451±0. 000012 with 
FWHM = 3.4±0.9 and 3.1±1.6 km s"^ respectively. These 
values translate to upper limits on kinetic temperatures of 
252 and 210 K respectively. The reduced in fliis case is 
0.87. This, as well as the residuals seen in Fig. [1] suggest 
that the two component model is probably over-fitting the 
data. 

The H I column density of an optically thin cloud cov- 
ering a fraction fc of a background radio source is related to 
the 21-cm optical depth t{v) in a velocity interval {v, v+dv) 
and to the spin temperature (Ts) by 

iV(H I) = 1.835 X 10^* / t{v) dv cm'^ (1) 

The total optical depth is / T{v)dv = 2.08±0.17 
km s~^. Thus the total H I column density is, 
A^(H l) = (3.82±0.31)xl0^*(rs//c)cm-^ From our optical 
spectrum we measure log iV(H l) = 21.36±0.10. For /c = 1, 
this yields Ts — 600t.^^g K. Interestingly this is very close 
to the upper limit we get from the width of the 21-cm ab- 
sorption line for the single gaussian fit. 

There is no milli-arcsecond scale resolution radio map 
available for the J 1337-1-3152 field. The radio source is 
compact in our GMRT images at 340 MHz with a resolu- 
tion of 10" x 8" and has a flux density of 69.3 mjy. In the 
VLA FIRST 1.4 GHz image (of resolution 5" x 5") the radio 
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Figure 2. Fit to the Lyman series lines of the z^ba = 3.17448 {v = km s ^) and 3.16768 {v = —488 km s ^) DLA systems. There is a 
strong Lyman hmit system at z^bs — 

3.1408 (v = -2420 km s'^). Some of the Lyman series absorption lines of the ^abs — 

3.17448 and 

3.1408 system happen to be blended. The vertical dashed lines at v ~ —1500 km s~^mark the locations of a C IV system along this line 
of sight (see also Fig. [Sj . 



source is represented by a single component with flux density 
83 mjy and a deconvolved size of 0.75" x 0.12". It is compact 
as well in our GMRT 610-MHz image with similar resolu- 
tion with a flux density of ~102mJy. The image of highest 
resolution for this source is from the JVAS/CLASS survey. 
In this VLA A-array 8 GHz image of resolution 0.26" x 0.22" 
the source is still compact with a flux density of 25.8 mJy. In 
the GB6 catalog at 5 GHz this radio source has a flux den- 
sity of 29mJy, very similar to the CLASS/JVAS 8 GHz im- 
age. Therefore, from the arcsecond and sub-arcsecond scale 
images available, the source appears to be core dominated 
with an indication of the low-frequency turnover. We have 
no reason to believe the covering factor is not close to unity. 

The spin temperature measured here is in between the 
expected kinetic temperatures o f cold (< 300 K) an d warm 
(5000-8000 K) neutral media (|Wolflre et al.l 120031 '). Thus, 
what we measure can be interpreted as the harmonic mean 



temperature of the multiphase medium along the line of 
sight. Interestingly the spin temperature we measure here 
is consistent with the mean spin te mperature measured in 
very low redshift QSO-galaxy pairs l|Carilli fc van GorkomI 
Il992h . 



4.2 Molecular Hydrogen 

Weak Lyman and Werner band H2 absorption lines from 
different rotational levels of the vibrational ground state are 
detected in system- 1. The spectrum at the positions of dif- 
ferent H2 absorption lines together with single component 
Voigt profile fits to these lines are shown in Fig. |3l All the 
lines except H2W3R1 and H2L5R0 were used during the fits 
for which z and b were assumed to be the same for all J 
levels. The reduced for the best fit is 1.4. Corresponding 
parameters are summarised in Table [1] 
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Table 2. Results of Voigt profile fits to the metal lines in system- 1 
at Zabs= 3.17447 
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Figure 3. Voigt profile fits to H2 Lyman and Werner band ab- 
sorption lines. The vertical dotted and dashed lines mark the 
locations of the 21-cm and H2 absorption lines, respectively. The 
apparent shift between the two lines is ~2.7 km s~^. The zero of 
the velocity scale is defined with respect to the redshift of the H2 
fines, Zabs= 3.174441. 



Table 1. Results of single component fits to H2 lines in system- 1 



Level z 


h 


log Ar(H2,J) 




(km s-i) 


(cm^^) 


J=0 3.174441(5) 


4.6±0.4 


12.88±0.16 


J=l 




13.85±0.04 


J=2 




13.65±0.04 


J=3 




< 13.47 









For J = 0, most of the lines are blended in the Lyman-Q 
forest and the quoted column density is constrained by the 
strength of the H2W1R0 absorption. The large error is due 
to the fact that this Une is blended with H2W1R1. Similarly, 
in the case of J = 3, most of the strongest transitions are 
blended with other lines in the Lyman-Q forest. We use the 
expected position of the H2W3R3 line to derive an upper 
limit on Ar(H2,J=3). 

From the total H2 column density and the above mea- 
sured A'^(H l) we derive the molecular fraction in the cloud: 
/h2 ~ 2iV(H2)/iV(H i) = 1.0 X 10"^ w hich is the lowest 
value measured in a DLA (see Table 1 of iNoterdaeme et al.l 
([2008a)). The low could be the result of either or both 
low dust content and excess UV flux from the QSO. 

The J = 1 excitation temperature is usually used as an 
estimator of the kinetic temperature. However, it is known 
that this estimator is robust only w hen the tot al H2 column 
density is in excess of 10 cm" (|Rov et al.i r2006'). which 
is far from being the case here. Indeed, the column density 
ratio A'^(H2,J=l)/Af(H2,J=0) measured here is higher than, 
or close to 9, the maximum value allowed by the Boltzmann 
distribution. Note that this has already been observed in the 
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log N 


h 
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(cm-2) 


(km s-l) 





N I 


3.174025 


13.78±0.03 


11. 3± 0.6 


-0.95±0.05 


q; n 

01 II 




i4.DU±U.Uy 




-t-U.26±U.lU 


S II 




13.98±0.04 






Fe II 




13.90±0.02 




-0.38±0.09 


C II* 


3.174387 


13.34±0.04 


8.3± 0.4 




C I 




< 12.30 






N I 




14.68±0.02 




-0.87±0.03 


I 




16.25±0.10 




-0.08±0.10 


Mg II 




15.13±0.06 




-0.05±0.06 


Si II 




15.17±0.04 




-l-0.00±0.04 


S II 




14.80±0.02 






Ar I 




13.62±0.06 




-0.38±0.10 


Cr II 




13.00±0.17 




-0.29±0.17 


Fe II 




14.53±0.03 




— 0.57±0.08 


Ni II 




13.41±0.04 




-0.44±0.04 


Zn II 




12.17±0.25 




-0.10±0.25 


C II* 


3.174473 


13.27±0.07 


2.5± 0.6 




C I 




< 12.30 
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-1.38±0.26 


I 




16.60±0.13 




+0.32±0.16 


Mg II 




15.03±0.06 




-0.10±0.22 


Si II 




15.11±0.11 




-f0.00±0.24 


S II 




14.75±0.21 






Ar I 




13.70±0.08 




-0.25±0.23 


Cr II 




12.89±0.22 




-0.35±0.30 


Fe II 




14.61±0.07 




-0.44±0.23 


Ni II 




13.26±0.05 




-0.54±0.22 


Zn II 




<12.30 







" Solar photospheric abundances are from lAsplund et al. 1 1 I2OO9II . 



low-Ar(H2) Zabs = 3 . 0248 system toward Q 0347-383 (see 
iLedoux et all l|2003D , ISrianand et al.l (|2005l )). 

From Fig.[3]it is apparent that, while the 21-cm absorp- 
tion is well within the velocity range of the H2 absorption, 
the centroid of the H2 absorption lines is shifted by 2.7±1.0 
km s~'^with respect to the centroid of the 21-cm absorption 
line for the single component fit. If we use the two compo- 
nent gaussian fit to the 21-cm line then the redshift of the 
weak 21-cm absorption component is consistent with that of 
the II2 absorptions within measurement uncertainties. The 
basic observation remains, however, that the optical depth 
ratio between H2 and 21-cm absorption varies inside the ab- 
sorbing cloud. 



4.3 Heavy element abundances and dust content 

We performed simultaneous Voigt profile fits to metal ab- 
sorption line features from neutral and singly ionized species 
(see Fig. |4)). The best fit is obtained with three components 
at Zabs = 3.174025, 3.174387 and 3.174473. First two com- 
ponents have relative velocities compared to the third one 
of A« — —32 and —6 km s~^ respectively. Fit parameters 
are summarized in Table [5] where the last column gives rel- 
ative abundances with respect to Sulfur in each component 
without applying any ionization correction. 
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Figure 4. Voigt profile fits to a selected set of metal absorption lines from system-1 at 2:abs= 3.17447. Vertical tick marks indicate the 
locations of the three velocity components. 



In Fig. m we also plot the 21-cm absorption profile for 
comparison. It is clear from this figure that the strongest 
and narrow metal line component at Zabs = 3.174473(6) is 
well aligned with the 21-cm absorption. It is also clear from 
Table [2] that there is no strong component to component 
variation in the depletion pattern among the three compo- 
nents which is consistent wit h what is observed in typical 
DLAs jRodn'guez et al.ll2006l ). 

Based on the robust column densities of Si ll, S ll and 
Fe 11 we find the main component at z^bs = 3.174473 to 
contain 40 to 45% of the total column density detected in 
this DLA. As there is no reason to believe the metallicity in 
this component is less than in the other two components, it 
is most likely that only 40 to 45% of the total H I column 
density is associated with this component. In that case we 
derive Ts ~ 270I72" ^ f^o'" ^'^is component. This brings the 
spin temperature very close to what is expected in a cold 
neutral medium (CNM). It is also interesting to note that 



there is no indication of a distinct narrow metal component 
coinciding with the II2 position supporting the fact that the 
cloud is not strictly homogeneous and that the H2 and 21-cm 
narrow components are consequences of physical conditions 
in the cloud changing on small scale. 

Adding column densities from the three components, we 
derive the average metallicities in the DLA given in Table [S] 
The a element metallicity is [Si/H] ~ [S/H] = -1.45 ±0.22. 
The oxygen abundance measured from the unblended lines 
matches well (i.e [0/S] = 0.05±0.29) with that of other a 
elements although the errors are large. 

The nitrogen abundance measurement is robust and 
we find [N/a] ~ —1.0. This is co nsistent with what is 
seen in low metalli city DLAs (see iPetitiean et al. I I2OO8I : 
iPettini et all l2008bl ) and very different o f what is usu- 
ally measured in gas associated with QSOs (| Pet it lean et al.l 
1 19941 : iHamann fc FerlandllTool ') . 

From Table [3l it can be seen that Ar seems slightly 
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Figure 5. Velocity plot of metal absorption lines of fiighly ionized species. The velocity scale is defined with respect to the DLA at 
2abs = 3.17447. Dashed vertical lines mark the locations of C IV components. The system at v ~— 500 km s~^ corresponds to the DLA 
a-t 2abs = 3.16768. The vertical dashed lines around v ~ —2420 km s~^are the C IV components of a LLS at Zabs= 3.1408. The C IV 
components marked at ~ —1500 km s~^are from an intervening high column density Lyman-o system. 



(~0.4 dex) underabundant compared to other a elements. 
Such deficiency has been noted already in the Galactic ISM 
and i n high redshift DLAs ijVladilo et al.ll2003l : [Ledoux et"all 
l2006f ). This has been attributed t o the lar^e photoioniz ation 
cross-section of Ar I above 15 eV l|Sofia fc Jenkinj|l998h . We 
can use this argument to constrain the ionization parameter 
in the absorbing gas (see below). Since the [Ar/Si] ratio is 
about the same in the individual components (see Table [2]), 
we can infer that the effect of ionization is nearly identical 
in the different components. 

The Zn ll transitions are weak and the corresponding 
column densities are difficult to measure. The [Zn/Cr] metal- 
licity ratio is close to solar while Fe, Cr and Ni are slightly 
underabundant compared to S and Si. If we attribute this 
to dust depletion, then we can infer that a maximum of 60% 
of Fe has gone into dust. In turn, the very low abundance 
of H2 may indicate that the gas is dust free with a slight 
enhancement of a element abundances. The chemical prop- 
erties of this system are very similar to the Zabs= 2.402 DLA 



towards HE 0027—1836 , that also shows low me tallicity and 
low molecular fraction l|Noterdaeme et aLll2007l ). 

From Fig. [5] it can be seen that C iv and Si iv absorp- 
tion lines are very strong in this DLA. N v is absent even 
though the absorption redshift is close to the QSO emission 
redshift. The absorption lines are seen at the expected posi- 
tions of C III, N III, Si III and O vi transitions. However, as 
these transitions are redshifted within the Lyman-a forest 
it is very difficult to conffrm them. This indicates however 
that a wide range of ionization states coexist in this system. 
Absence of N v and the possible presence of O vi is con- 
sisten t with what is typically seen in DLAs (see IFox et all 
l2007f) . 



4.4 Physical conditions in the absorbing gas 

Here, we discuss the physical conditions in the gas asso- 
ciated with syst em- 1 using simple p hotoionization models 
using CLOUDY (|Ferland et al.lll998f V The aim is not to de- 
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Table 3. Average metallicities in system-1 at Zabs = 3.17447 



Species 


log N 


log Z/Zq" 









H I 


21.36±0.10 




C I 


< 12.60 




C II* 


13.61±0.08 




N I 


14.83±0.06 


-2.44±0.12 


O I 


16.76±0.16 


-1.40±0.19 


Mg II 


15.59±0.08 


-1.35±0.13 


Si II 


15.50±0.15 


-1.42±0.18 


S II 


15.11±0.20 


-1.45±0.22 


Ar I 


13.96±0.10 


-1.80±0.14 


Cr II 


13.25±0.28 


-1.79±0.30 


Fe II 


14.91±0.08 


-1.95±0.13 


Ni II 


13.64±0.15 


-1.97±0.18 


Zn II 


12.26±0.26 


-1.77±0.28 



°' solar photospheric abundances are from I Asplund et al. I 1I2OO9I') . 
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log Ionization paramotcr 



scribe in detail the system but rather to estimate the ion- 
ization parameter using A'^(Ar l), the particle density using 
A'^(C II*) and to constrain the distance of the absorbing gas 
from the QSO. 

We assume the absorbing cloud to be a plane paral- 
lel slab of constant density, the ionizing radiation to be 
that of a typical QSO (i.e iMathews fc Ferlandl l|l987l ')'). 
log Af(H i) = 21.3 and Z = -1.4 Zq. The observed dust 
to metal ratio in the system is ~ 60% of the Galactic value. 
Together with the measured metallicity we find the dust to 
gas ratio to be roughly 2% of wh at is measured in the cold 
phase of the neutral galactic ISM l|Weltv et al.lll999l 'l . We use 
this value in our model and assume the dust composition to 
be similar to that in the galactic ISM. 



Constraints on the ionization parameter 

In Fig. [B] we plot the Ar l/Ar ratio as a function of the ioniza- 
tion parameter, U. Assuming solar [Ar/S] (and/or [Ar/Si]) 
abundance ratio, we derive from the Ai'(Ar l)/N{S 11) 
(and/or N{At l)/iV(Si 11)) ratio given in Table [3] that the 
Ar^/Ar ratio is larger than 0.4. From Fig. [S] we conclude 
U < 10~^'^. Note that as the average Ar^/Ar ratios are 
similar in different components, this result will hold for in- 
dividual components as well. 

In principle we can derive a lower limit on the ioniza- 
tion parameter from the upper limit on A^(C l). As the C 11 
absorption lines available in our spectrum are heavily sat- 
urated we have to a ssume a value o f the [ C/Si] abundance 
ratio. From Fig. 3 of iPetitiean et all (|200^ we note that Si 
and O abundances follow each other very well over a wide 
range of metallicities. However, carbon can be underabun- 
dant (by ~ 0.5 dex) compared to oxygen for [O/H] ~ —1.4 
as observed in this system. Using [C/Si] — —0.5, we derive 
log 7V(C 11) = 15.96±0.15 (i.e log N{Si n)+8.52-7.56-0.50). 
This, together with the 3a upper limit on A''(C l), gives 
N{C i)/N{C II) < 10"^ Fig. Ogives the Ai'(C i)/iV(C 11) 
ratio through the cloud model versus ionization parameter. 
The above condition implies log U > —3.3. We therefore find 
—3.3 <log U < —3.1 and we will take log U — —3.2 as a 
typical value for further discussions. 



Figure 6. Fraction of neutral Ar as a function of the ionization 
parameter. 
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Figure 7. The predicted A''(C l)/A'"(C 11) ratio as a function of 
ionization parameter. 



Fine- structure lines and particle density 

The only fine-structure absorption line seen in this system 
is C n*A1335. Using the same assumptions as above, we 
can derive log iV(C ll*)/Af(C 11) = -2.45 ± 0.16 from the 
observations. Interestingly this is consistent with what is 
expected in the case of CNM gas with low metallicity (See 
fig. 15 of lSrianand et al.ll2005l ). 

In the absence of C l and associated fine-structure lines 
it is not easy to derive nu as the C 11* column density 
depends on the ionization fraction of the gas and the ki- 
netic temperature (see, for example, the discussion in sec- 
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Figure 8. The slanted solid line gives the log Af(C n*)/N(C ll) 
ratio as a function hydrogen density. The observed mean and Icr 
ranges for N{C ll*)/N{C ll) are given by horizontal solid line and 
dashed lines, respectively. The intersections of these lines with 
the model curve give the allowed density range. 

tion 6 of ISrianand et al]|2005l ). Therefore, to interpret the 
N{C u'')/N{C ll) ratio, we will rely on a photoionization 
model. We run a grid of Cloudy models keeping constant the 
ionization parameter, log U = —3.2, as constrained above, 
T — 600 K, as derived from the 21-cm spin temperature 
measurements, and varying the particle density in the range 
0.1 to 10 cm"^. 

The model predictions on A'^(C n*)/7V(C ll) are given in 
Fig. [8] as a function of the particle density. We derive 1.3 < 
riH (cm~'^) < 2.6. Thus a value of 2 cm"'' is representative 
of the particle density in this system when we assume log 
U — —3.3. Cloudy also predicts the column densities of the 
O I and Si ll fine-structure lines and the predicted values 
are consistent with the non-detections of these species in 
our data. 

As an aside, we notice that while most of the constraints 
based on atomic species are satisfied by the model presented 
here, the latter predicts high values of 7V(H2) (remember 
we use a dust to gas ratio equal to 2% of the Galactic ISM 
value). However, models with very low dust content (i.e dust 
to gas ratio <0.1% of the Galactic ISM) predict an amount 
of H2 consistent with our observations. Thus it may be pos- 
sible that the observed lower abundance of Fe co-production 
elements does not refiect dust depletion but real nucleosyn- 
thetic underabundance. 

4.4-3 Constraints on the distance to the quasar 

From the estimates of the ionization parameter and parti- 
cle density, we can constrain the distance of the absorbing 
cloud from the QSO. The Lyman limit of the absorber in 
our observed frame is at A ~3807A. From the fiux calibrated 
SDSS spectrum we find fu — 70 /iJy at this wavelength. For 
a flat cosmological model with JIa = 0.73, fim ~ 0.27 and 



h = 0.7 we find the luminosity distance of the QSO to be 
2.8x10^* cm. This gives the luminosity at the Lyman limit, 
L„{LLS) = 1.7x10^° ergss"^ Hz"^ and the rate of hydrogen 
ionizing photons, (5(H) — 2.6x10^^ photons s~^, when we 
assume a flat spectrum. The ionization parameter is usually 
deflned as. 



Here, r^ is the distance of the cloud from the QSO and n is 
the particle density in the absorbing gas. Using n ~ 2 cm""' 
and U = 5xl0~* we get rc ~ 270 kpc. This also means that, 
if in situ star formation controls the ionization state of the 
gas then the DLA has to be at a distance much larger than 
270 kpc from the QSO. 

This clearly shows that the absorbing gas is not as- 
sociated with the QSO host galaxy. The gas could be as- 
sociated with a galaxy in the group or cluster in which 
the QSO resides. In this sense this system is very similar 
to that of the 2ab3= 2.8112 DLA towards PKS 05 28-250 
(ISrianand fc Petitiean|[l998l : iMoller fc Warrenlll998l ). 

4.5 Constraining fundamental constants 

As the energy of the hyperfine H I 21-cm transition is propor- 
tional to the combination of three fundamental constants, 
X = a^Gp/fi, high resolution optical and 21-cm spectra can 
be used together to pr obe the combined cosm ological varia- 
tion of these constants (|Tubbs fc Wolfelll980l). For the 2abs = 
1.7764 system towards Q 1337-1-170. ICowie fc Songaili 
(jl995l ) have obtained Ax/x = (2uv-22i)/(l + 22i) = (0.7± 
1.1) x 10~^. Here, zuv and 221 are the re dshifts of the UV 
l ines and 21-cm absorption, respectively. iTzanavaris et all 
(2005) obtained Ax/x = (0.63 ± 0.99) x lO"'^ using this 
technique and a sam ple of 8 published systems. Similarly, 
iKanekar et al] (|2006l ) derived Ax/x = (-0.15±0.62) x 10"^ 
at Zabs= 2.34786 towards PKS B0438-438. 

The measured wavelengths of H2 Lyman and Werner 
band tra nsitions can be used to probe the variation of /i 
( V arshalovich fc Levshakovlll993h . The presence of H2 to- 
gether with 21-cm absorptions can, in principle, allow us 
to place constraints on a and fi in the same system. In 
this regard, the system studied here is unique. However, 
as the H2 column density is weak, only Werner band ab- 
sorption lines are seen. Therefore and unfortunately, the 
range of sensitivity coefficients of the detected lines is 
too narrow to provide stringent constraints. Using a lin- 
ear regression analysis, we derive Afi/fi < 4.0 x 10"* and 
ZH2 = 3.174442±0.000014. The upper limit on A^i/^i mea- 
sured here is roughly an order of magnitude larger than 
the typical estimates ba s ed on H2 lines publis h ed till now 
[see'lvanchik et al.' (2005') 'Reinhol d et~al] (|2006l '). lKing et all 
(2008) and Thompson et al. (200!|)]. 

As discussed before, the single component fit to 
the 21-cm absorption gives 221 = 3.174480±0.000002. 
The strongest metal line component is seen at 
zuv = 3. 174473±0. 000006. From these two redshifts, 
we derive Ax/x = -(1.7 ± 1.5) x 10"®. We find the 
wavelength calibration accuracy is < 4 mA duri ng the 
pipeline data reduction (see also IChand et al.l l2006f ) . This 
translates to an error of 4xl0~® (or a velocity uncer- 
tainty of 0.3 km s~^) in the redshift measurement at the 
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mean observed wavelength (~ 4500A) of the lines used 
to derive zuv- If we include this error also then we get 
Ai/a: = -(1.7 ± 1.7) x lO"''. 

If we assume no variation in a and Gp (respectively 
^ and Gp) then this limit translates to A/x//i = — (1.7± 
1.7) X lO"'^ (respectively Aa/a = -(0.85 ± 0.85) x lO"**). 
These constraints alone are more stringent than most of the 
published values in the literature. 

However, we like to point out that our UVES spectrum 
was taken as part of our ongoing survey for molecules in 
DLAs where ThAr calibration frames were not taken im- 
mediately after the science observations. This may induce 
an additional systematic error of the same order as our 
measurement errors. In addition, the velocity offset we no- 
ticed between H2 and 21-cm absorption could reflect inho- 
mogeneities in the gas which could contribute to additional 
systematic error. This kind of systematics can only be alle- 
viated through randomization using a large and well defined 
sample that we want to build in the future. Note that the 
velocity shift could be related to covering factor issues if the 
radio source shows structures at the milli-arcsec scale. Thus, 
as for any other system, high resolution VLBA observations 
are needed to support the analysis towards constraining vari- 
ations in X. 



o 
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Figure 9. Multicomponent Voigt profile fits to absorption fea- 
tures in the 2abs= 3.16768 system. Vertical dotted lines mark the 
locations of individual components. 



5 ANALYSIS OF SYSTEM-2 (z^bs = 3.16768) 

Molecular hydrogen and C ll* absorption lines are not de- 
tected in this DLA with log iV(H l) = 20.41±0.15. The 2a 
upper limit on the molecular fraction is fn^ < 2xl0~^. In 
addition, as our GMRT observations were performed prior 
to our UVES observations, we did not cover the expected 
frequency range for the redshifted 21-cm line. However, this 
system is very interesting as the metal lines from strong res- 
onance transitions are unsaturated and D I absorption can 
be fitted. This makes it possible for us to understand the 
chemical evolution and deuterium astration factor in this 
system which is at an early stage of chemical evolution. 



5.1 Heavy element abundances and dust 

Low ionization metal absorption lines detected in this sys- 
tem together with the best fitted Voigt profiles are shown 
in Fig. [51 The strongest transitions of C ll. Si ll and O I 
are not strongly saturated and can be used to derive col- 
umn densities. The absorption profiles are well fitted with 
four components at « = —7.7, 0, 12.4 and 17.8 km s~^ with 
respect to Zshs= 3.16768. 

Table[4]summarises the total column densities of neutral 
and low ionization species and the corresponding metallic- 
ities measured in this system assuming no ionization cor- 
rection. The oxygen metallicity is one of the lowest know n 
among DLA s and s ub-DLAs (see iPetitiean et al.l l|2008f ): 
iPettini et all l|2008bl )). We remind the reader that the low 
oxygen metallicity cannot be due to ionization effects be- 
cause O I is known to be coupled tightly with H I through 
charge exchange reactions. Actually the carbon, silicon, alu- 
minium and iron metallicities estimated from the singly ion- 
ized species are consistent with [0/H]. 

The system is of low-ionization. We detect the 
Si inA1206 fine and measure log 7V(Si ill) = 13.28±0.06, but 



Table 4. Total column density and average metallicity in the 
^abs = 3.16768 DLA 



Species 



log N 
(cm^^) 



z/z- 



H I 


20.41±0.15 




N I 


<12.80 


< -3.53 


O I 


14.43±0.09 


-2.64 ±0.17 


C II 


13.98±0.06 


-2.82 ± 0.16 


Si II 


13.24±0.05 


-2.68 ± 0.16 


Si III 


13.28±0.06 




Al II 


12.00±0.05 


-2.78 ± 0.16 


Fe II 


13.14±0.26 


-2.72 ± 0.30 



solar photospheric abundances are from Asplund et al. (2009) 



the Si IV absorption is very weak (see Fig. [5]). Therefore, in 
any case the silicon abundance is [Si/H] < —2.38 ± 0.13. 

The spectrum covers the expected position of C lllA987 
but this line is blended with intervening Lyman-a absorp- 
tions (see Fig. [5]). C iv absorption is also weak. The column 
density of C ll implies [C/0]> -0.18 ±0.18. This value con- 
tradicts the result of standa r d gal actic chemical evolution 
models (see lAkerman et al.l (|2004l ')') expecting the carbon 
metallicity to be much smaller than the oxygen metallic- 
ity (i.e [C/O] < -0.5) for [ 0/H] < -2. H o wever our result 
is consistent with result s by Pettini et al.l 1 2008bl ) and sup- 
ports the suggestion bv lAkerman et aL 1 l|2004l ) that carbon 
and oxygen may be produced in near-solar proportions in the 
earliest stages of galactic chemical evolut ion. Chemical evo- 
l ution models with Population III stars (jChieffi fc Limongil 
(i2003 )l or with fast rotating stars (|Chiappini et al.l (120061 )) 
produce near solar [C/O] at low metallicities. 

The spectrum does not show detectable absorption lines 
of N I. Hence, we only derive an upper limit for N(N l). 
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Figure 10. Voigt profile fits to the Lyman series lines in the 2abs= 3.16768 system. Ticks around ti ~ km s~^mark the locations of the 
metal line components. The Voigt profile parameters giving good fits to the damping wings of Lyman-/3, the red wings of Lyman-7 and S 
and the blue wings of Lyman-6, Lyman-7 and Lyman-8 transitions (profiles with long dashed lines) do not explain the excess absorption 
at the expected positions of the D I absorption lines (ticks around v = —82 km s~^). Fits to D I absorption lines are shown with dashed 
profiles. 



The expected -wavelength range for N nA1083 is affected 
by intervening Lyman-a absorption from the IGM. Using 
N l/O I -we derive [N/O]<-0.89. This is not stringent enough 
to draw conclusions regarding the contribution to the nitro- 
gen metaUicity from primary and secondary productions. 

Based on the [Fe/a] abundance ratio we can conclude 
that there is no dust depletion in the gas. The absence of 
H2 absorption could very well be attributed to the lack of 
dust in the system. The lack of a element enhancement also 
means that the star formation activity is not recent. In addi- 
tion, this also undermines the contribution of Population III 
stars to overall enrichment. 



5.2 Constraint on the D/H ratio 

While fitting the Lyman series absorption lines, we find a 
clear excess absorption in the blue wing of all the available 
Lyman transitions (see Fig. I10|l . The total H I column den- 
sity is well constrained by the damping wings of Lyman-/? 
and Lyman-7. The uncertainty in its measurement comes 
from the error propagated through the Voigt profile fits to 
the Lyman-/? absorption of System- 1 and from continuum 
placement. The lower panel of Fig. [10] also shows a typi- 
cal continuum used near the wavelength range covered by 
Lyman-/3 line. This wa s determined using the S PSS com- 
posite QSO spectrum l|Vanden Berk et al.l I2OOII ) with ap- 
propriate scaling to fit the line free regions. We constrain 
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Figure 11. Top panel: D/H values measured in absorption sys- 
tems as a function of A'^(H l). Bottom panel: D/H as a function 
of metallicity. In both panels the open circle is the measurement 
presented in this paper. Filled circles are earlier measurements 
ba sed on N(T) l)/Af(H l ) lines from the literature (as summarised 
m llvanchik et al.ll2010l ). Stars are the measurements based on 
HD /2H2 in two DLAs p resented in iNoterdaeme et all ll2008bl ') 
and llvanchik et alj 1I2OIC1I) . 



the velocity dispersion by using the red wings of Lyman-7 
and 5 and using the blue wings of Lyman-6, Lyman-7 and 
Lyman-8. We notice that the resulting best fitting Voigt pro- 
files (dotted profiles in Fig. llOp under-predict the blue wing 
absorption at ~— 82 km s^^for Lyman-7, ^ s-nd ^- Consis- 
tent absorption dips are present at similar velocity in the 
blue sides of higher Lyman series lines. We do not find any 
metal line absorption at this velocity (see Fig. It is also 
clear from Fig. [2] that the excess can not be explained by 
blending with Lyman series lines from the other two strong 
H I absorption systems seen in the immediate vicinity. Thus 
it is most likely that the excess absorption in the blue wings 
of H I lines is due to D i. 

We simultaneously fit the D I Lyman series lines with 
three components at the redshifts of the three strongest 
O I components that contain ~ 96% of the O I column 



density. We varied A'^ and b for the D I lines to get the 
best fit. It is to be noted that D 1—6 line is unsaturated 
and the total A^(D l) measurement is insensitive to the ex- 
act component structure. The results of the Voigt pro- 
file fits are shown in Fig. 1101 The best fit model gives 
A^(D i) = (3.0 ± 0.2) X 10^^ cm-2 leading to log (D/H) = 
log (D i/H l) = -4.93 ±0.15. This is a factor of 2 lower than 
the primordial value, log [D/H]p = —4.59 ± 0.02, derived 
from five-year d ata of the Wilkinson M icrowave Anisotropy 
Probe (WMAPI llKomatsu et ahlboogl ). 

For the metallicity estimated in this system, it is be- 
lieved that [D/H] will be close to [D/H]p. However, the 
present system seems to have lost an appreciable fraction 
of deuterium during its evolution. In Fig. [11] we plot D/H 
measured in this system together with those from the liter- 
ature as a function of 7V(H l) and [O/H]. The li t erature 
data includes the systems listed b y iPettini et al.l (|2008a ) 
together with the measurement by ICrighton et al] (2004) 
in a Lyman limit absorption system at 2: = 3.256 to- 
wards PKS 1937— 1009 and the measurement in the ZabB = 
3.02 DLA towards 0347-3819 (|Levshakov et all 120021 '). For 
this system, we use the [O/H] measurement given in 
iPetitiean et al.l (|2008l ). We a lso plot the D/H measurem ents 
based on the HD/H2 ratio bv lNoterdaeme et al.l l|2008b! ) and 
llvanchik et al.l (|2010l ). In both panels of Fig. [11] the horizon- 
tal lines give the we ighted mean and Icr range reported by 
iPettini et all l|2008al ) using their sample of 7 measurements. 

The latter authors identify t wo outliers, the ZabB = 
2.07623 DLA towards Q 2206-199 l|Pettini fc Bowenll200ll ) 
and the Zaba = 2.50357 partial Lym an limit system towards 
Q 1009-1-299 (|Burles fc TytlerlllOTSj l. that may be responsi- 
ble for the scatter in their sample. Interestingly, the value 
measured here is consistent with log (D/H) = —4.78 ± 0.09 
measured in the Zabs= 2.07623 metal poor DLA towards 
Q 2206—199. Therefore, addition of our measurement to the 
existing sample increases the scatter (see Fig, illf) . It is also 
interesting to note that a trend seems to exist for higher D/H 
values at higher oxygen metallicity which may be difficult 
to understand. 

As the absorption lines of H I and D I fall inside the 
Lyman-Q forest one can always question the accuracy of 
column density measurements. As we use damping wings, 
the A'^(H l) measurement is robust as is our lower limit on the 
astration factor. Indeed, if the deuterium absorption features 
are due to intervening Lyman-a clouds, the astration factor 
must be even larger. 

There are recent indications that infalling gas may cause 
astration factor to be different fr om that predicted by chem - 
ical evolution closed box models. INoterdaeme et ahl (|2008bl ) 
measured a relatively low astration factor of deuterium cou- 
pled with a high metal enrichment in t he 2a ba= 2.418 DLA 
towards Q 1439-1-1117. llvanchik et al.l (|2010l 'l have reported 
a [D/H] value close to [D/H]p for the «abs= 2.338 DLA to- 
wards Q 1232-1-0815 based on H2 and HD column densitie s. 
This system has [O/H] = -1.43±0.08 ('Srian and et"al]|2000l '). 
Thus it appears that the astration factor of deuterium varies 
largely from one system to the other at a given metallicity. 
Therefore, it may not be a good idea to remove the so-called 
outliers while computing the weighted mean value of [D/H] 
measured in absorption line systems. 
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6 CONCLUSION 

We have presented detailed analysis of two very closely 
spaced DLAs (at Zabs= 3.17447 and 3.16768) within a ve- 
locity separation of ~500 km s^^along the line of sight to 
SDSS J1337+3152. We have reported the detection of 21- 
cm and molecular hydrogen absorption lines in the damped 
Lyman-Q system at 2abs= 3.17447. This is the second time 
that both species are detected in the same system. This is a 
unique combination to constrain at the same time the time 
variation of the fine structure constant, a, and the proton- 
to-electron mass ratio, /x. We note however a velocity shift 
of ^^2.7 km/s between the H2 and 21-cm absorption, reveal- 
ing small scale inhomogeneities in the cloud that could wash 
out detection of a variation in x. Such systematics can only 
be alleviated with a large sample of such systems. However, 
the redshift of the strongest metal line component matches 
well with the redshifted 21-cm absorption line. Using this, 
we derive a stringent limit on the variation of j: = a^gp/fi, 
Ax/x = -(1.7 ± 1.7) X 10"", which is one of the best limits 
ever obtained. While this demonstrates the advantages of 
combining GMRT and UVES data, we emphasize the need 
for a large sample of such measurements. 

Using A^(H l) measured from Voigt profile fitting of the 
damped Lyman-a line in the UVES spectrum together with 
the 21-cm optical depth measured from our GMRT spec- 
trum, we obtain Ts/fc = GOOtlH K. This, together with 
the velocity shift noted above, is consistent with the absorb- 
ing gas being a mixture of different phases. We derive mean 
physical conditions in the cloud, U ~ 10~^'^, nn ~ 2 cm~"^, 
Tk ~ 600 K. The metallicity is of the order of [S/H] ~ -1.5 
and the distance to the quasar larger than 270 kpc. We con- 
clude that the gas must be associated with a galaxy in the 
group or cluster in which the QSO resides. 

The gas responsible for the Zabs= 3.16768 DLA has very 
low metallicity, [O/H] = -2.64 ± 0.17. This allows us to 
measure the column density of C 11 accurately. We mea- 
sure [C/O] = -0.18±0.18. This means that [C/0] reaches 
solar value even at very low metallicities. We measure the 
deuterium abundance in this DLA to be log (D/H) — 
log (D l/H l) = -4.93 ± 0.15. This is a factor of 2 lower 
than the primordial value, [D/H]p — —4.59 ± 0.02, derived 
from five-year dat a of the Wilkinson M icrowave Anisotropy 
Probe (WMAP) dKomatsu et al.ll2009il . We therefore con- 
clude that astration factors can vary significantly even at 
low metallicity. 
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